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Preface 


The majority of all stars are members of a binary system. The evolution of 
such binary stars and their subsequent production of pairs of compact objects 
in tight orbits, such as double neutron stars and double black holes, play a 
central role in modern astrophysics, Binary evolution leads to the formation 
of different types of violent cosmic events such as novae, supernova explosions, 
gamma-ray bursts, mass transfer and accretion processes in X-ray binaries, and 
the formation of exotic radio millisecond pulsars. In some cases, the binary 
systems terminate as spectacular collisions between neutron stars and/or black 
holes. These collisions lead to powerful emission of gravitational waves, as 
detected by LIGO since 2015. The coming decade is expected to reveal a large 
number of discoveries of binary compact systems, as well as their progenitors and 
merger remnants, from major instruments such as the radio Square-Kilometre 
Array; the gravitational wave observatories LIGO-Virgo-KAGRA-IndIGO and 
LISA; the astrometric space observatory Gaia; the James Webb Space Telescope; 
and the X-ray space observatories eXTP, STROBE-X, and Athena. In this 
light, it is important to have a modern textbook on the physics of binary stars 
evolution, from ordinary stars to X-ray binaries and gravitational wave sources. 
The scope of this book is that the reader (student or educated expert) will 
learn the physics of binary interactions, from stellar birth to compact objects, 
and relate this knowledge to the latest observations. The reader will learn about 
stellar structure and evolution, and detailed binary interactions covering a broad 
range of phenomena, including mass transfer and orbital evolution, formation 
and accretion onto compact objects (white dwarfs, neutron stars and black 
holes), and their observational properties. Exercises are provided throughout 
the book. 


It has been a privilege for us to write this book, and we hope you enjoy it. 


Thomas Tauris & Ed van den Heuvel — June 2022 
(Aalborg / Aarhus / Bonn) (Amsterdam) 


Figure 1: The authors (Ed van den Heuvel and Thomas Tauris) enjoying a 
relaxed evening in Thomas’ house in Djursland, Denmark, after a long day’s 
work on the book in February 2020. [Photo by Birgitte Tauris.] 


Chapter 1 


Introduction: 
The Role of Binary Star 
Evolution in Astrophysics 


Many key astrophysical objects and phenomena are related to the evolution 
of binary stars. This holds, for example, for the formation of the brightest 
X-ray sources in the sky and the formation of double neutron stars (NSs) and 
black holes (BHs), the mergers of which produce the strongest bursts of energy 
anywhere in the observable Universe, measurable on Earth with gravitational 
wave (GW) detectord"| Binary interactions also play a role for the origin of 
most supernovae (SNe), all nova explosions and short gamma-ray bursts, and 
the formation of millisecond radio pulsars and a large variety of stars with 
peculiar chemical abundances, such as barium stars and carbon-enriched metal 
poor stars. 

It has been long realized, from the fact that most stars are members of 
binary systems 2003), that binary evolution 
must play a key role in stellar evolution (e.g.|van den Heuvel| (1994). This early 
and important awareness has even been strengthened further by the findings of 
e.g. (2011); (2012), that practically all massive stars are 
found in binaries with orbits such that at some stage in their evolution, the far 
majority of these stars will interact with each other. This implies that binary 
interactions dominate the evolution of massive stars. 

The first ideas about the evolution of binary systems originated in the 1950s. 
They were largely inspired by the surprising characteristics of Algol-type eclips- 
ing binary systems (see also Section 2.4). Here, Algol-type means physically 
similar to the Algol-system, which consists of an unevolved B8V main-sequence 
star with a mass of 3.2 Mo together with an evolved but less massive subgiant 


1The first ever detected BH merger event GW150914 (Abbott et al.| |2016) released an 


energy of ~3 Mo c? within a fraction of a second, thereby outshining all stars in the Universe 
for a brief moment. 


companion star of spectral type G8III, of mass 0.7 Mo. This situation, with 
the more evolved star having the smaller mass of the two, is just opposite to 
what one would expect on the grounds of stellar evolution, as stars of larger 
mass are expected to have shorter lives than stars of smaller mass do. In a 
binary — where both stars were born at the same time — one would therefore 
expect the star of larger mass to be in a more advanced stage of evolution at 
any time than that of its companion of smaller mass. 

This is what is called the Algol paradoz. was the first to 
realize that this paradoxical situation can be explained if one assumes that 
large-scale mass transfer can take place during the evolution of a binary system: 
Crawford hypothesized that the subgiant components in Algol-type binaries 
were originally the more massive components of these systems. As the more 
massive star evolved faster than its less massive companion did, it was the first 
one to exhaust the hydrogen fuel in its core and evolve into a giant star with 
a much expanded envelope. The presence of the close companion, however, 
prevented such an evolution: when the outer layers of the expanding (sub)giant 
came under the gravitational influence of the companion, they were captured 
by this smaller star (the accretor), causing the accretor to increase in mass at 
the expense of the (sub)giant (the donor). The (sub)giant transferred so much 
of its mass that it was finally able to restabilize its internal structure. At that 
moment, it had become the less massive of the two stars, and the originally less 
massive star became the most massive one of the pair. 

The first attempt to carry out a real calculation of this type of evolution 
with mass transfer was by [Morton| (1960). He demonstrated the correctness of 
Crawford’s conjecture, that mass transfer, once it begins, continues until the 
(sub)giant has become the less massive star of the system. In his calculations, 
however, Morton still assumed that the orbital period of the system does not 
change during the mass transfer. This is not correct because, if one assumes 
that the total mass of the system is conserved, one also expects the total angular 
momentum of the system to be conserved. The total angular momentum is, in 
good approximation, equal to the orbital angular momentum (as the rotational 
angular momentum of the two stars is usually much smaller than the orbital 
one). Conservation of the orbital angular momentum implies that, during the 
mass transfer the orbital period and separation change in a well-determined 
way, which will be explained in Chapter 4. The evolution of close binaries in 
this more realistic approach was first calculated, independently of one another, 
by [Paczyński] (1966), [Kippenhahn & Weigert (1967), and [Plavec] (1967). Their 
work was the foundation of all subsequent work on the evolution of close binary 
systems. Furthermore, the discovery of the first celestial X-ray source in 1962 
and of the X-ray binaries in 1971-1972 — which earned Riccardo Giacconi 
the 2002 Physics Nobel Prize — has given a great stimulus to the research in 
this field. The X-ray binaries consist of a normal star together with a compact 
object: a NS or a BH, which are the end states of evolution of massive stars. The 
X-rays from these binaries, which are observed with space-borne detectors, are 
generated by the accretion of matter onto the compact star, which is captured 
from the outer layers of the companion star. During the accretion this gas, 


falling inwards in the extremely strong gravitational field of the compact star, 
is heated by the release of gravitational potential energy to temperatures above 
10° K, causing it to emit X-rays. 

Without the occurrence of extensive mass transfer from the original primary 
star (the progenitor of the compact object) to the secondary star, most of these 
systems could not have survived the SN explosion of the primary star in which 
the NS or BH was formed (van den Heuvel & Heise) [1972 [Tatukov & Yangelson| 
(1973), because the probability of the post-SN orbit to remain bound depends 
on the relative mass loss from the system during the SN (see e.g. Section 4.3.10 
and Chapter 13). 

Since the early 1970s, the realization of the importance of accretion of matter 
onto a compact star (NS, BH or white dwarf [WD]) as an energy source in 
many types of binaries, ranging from X-ray binaries to cataclysmic variables 
(CVs) and symbiotic stars, has been a further important source of inspiration 
for new research on the structure and evolution of close binary systems. The 
discoveries since 1974 of binary radio pulsars, of which at least 20 are double 
NSs (DNSs; for a review, see|Tauris et al.||2017), have revealed many interesting 
properties, including the many relativistic effects that are measurable in them 
with unprecedented precision (c.g. [Taylor & Weisberg} {1989} {Taylor} |1992} [Kaspi] 
Bo). 


These discoveries have created a new and fundamental branch of relativistic 
binary star astrophysics, which among other things has produced the most ac- 
curate measurements of masses of any stellar objects so far (Chapter 14). The 
measured rate of orbital decay of the first-discovered DNS, PSR B1913+16, is 
in exact agreement with the decay rate predicted from the emission of GWs ac- 
cording to the general theory of relativity. The highly precise detection of this 
and other relativistic effects earned the discoverers of this binary pulsar, Russell 
Hulse and Joseph Taylor, the 1993 Physics Nobel Prize. Later-discovered DNSs, 
particularly the double pulsar system PSR J0737—3039, have further refined the 
up to five tests of Wat allowed by these systems (Section 14.7.1) to almost 
incredible precision [Wex] [2014] [Kramer et al.| (2021). 

The amazing “islets of the o a of a double BH, starting with 
GW150914 in September 2015 (Abbott et al.| [2016), and of a DNS, in Au- 
gust 2017 ( (Abbott et al. [2017b), have revealed the ultimate final destiny of a 
massive binary star system and demonstrated the production of strong bursts of 
GWs, observable on Earth. This earned the LIGO pioneers Rainer Weiss, Kip 
Thorne and Barry Barish the 2017 Physics Nobel Prize. But how do ordinary 
stars born in a binary system end up as two NSs or BHs, that finish as a final 
single BH remnant? To answer this question, we need to follow the binary sys- 
tem through a long chain of exotic binary interactions (Fig. 1.1), involving mass 
transfer between the stellar components, SNe, and relativistic effects. That the 
orbits of evolved massive stars in binaries will shrink to the very small sizes ob- 
served for the double compact objects, was predicted before these objects were 
discovered 1973). 

Radio pulsars are some of the most intriguing astrophysical objects, among 
which a sizable fraction of binaries is found with quite special characteristics that 
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Figure 1.1: Formation model of a close double neutron star (DNS) system as 
final product of the evolution of a massive close binary. The DNS system may 
eventually merge and leave a solitary BH remnant. The same model, scaled 
up to higher initial stellar masses, is one of the main scenarios to explain the 


formation of double BHs (see Chapters 10, 12, and 15). (After [Tauris et al.| 
2017). 
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Figure 1.2: Observed population of ~ 2400 radio and high-energy pulsars with 
measured values of both spin period (P) and period derivative (P). The vari- 
ety of pulsars has flourished immensely since their discovery in 1967 and con- 
tinues to be a major science driver in modern astrophysics. Binary pulsars 
(blue circles) dominate among the fast-spinning millisecond pulsars. CCOs 
are central compact objects of SN remnants, XDINSs are X-ray dim iso- 
lated NSs, and RRATs are rotating radio transients (see Chapter 14). Data 
taken from the ATNF Pulsar Catalogue in June 2021 (Manchester et al. 


https://www.atnf.csiro.au/research/pulsar/psrcat). 


give important information on binary evolution. To understand the continuously 
growing diversity of observed radio pulsars (see Fig. 1.2), it is necessary to 
link their properties to the stellar and binary evolution of their progenitors 
1901). 

The detection of close binary pulsars, and of the merger events of double BHs 
and NSs have further increased the interest in the evolution of binary systems. 
These discoveries demonstrated that, even though binaries may have undergone 
several stages of mass transfer during their evolution, where up to 90% of their 
original mass and > 95% of their original orbital angular momentum is lost from 


the system, and despite having experienced two SN explosions, the two stars 
might still survive as a (very close) binary system with two compact objects. 
The existence of the close DNSs and close double BHs has also demonstrated 
that a precise knowledge of the physics of binary evolution is of vital importance 
for understanding fundamental astrophysics as diverse as the generation of the 
strongest bursts of GW radiation, the production of gamma-ray bursts, and the 
synthesis of heavy r-process (or rapid neutron capture) elements in the Universe. 
The latter was predicted to be a ae by the merging of DNSs or NS+BH 


binaries (Lattimer & Schramm |1976) and expected to be observable as a so- 
called kilonova optical-infrared ora i (Metzger et al. | {2010} [Berger] [2014}. 


These predictions have been beautifully confirmed by the spectroscopic study 
of the optical-infrared transient that accompanied the DNS merger event 


GW170817 (Abbott et al.| |2017b}djc). The mergers of double compact object 


binaries with at least one NS also had been predicted to produce gamma-ray 


bursts (Paczynski||1986}|Eichler et al.|!1989), which was confirmed by the Fermi 


and INTEGRAL missions by the detection of a short gamma-ray burst, follow- 
ing within two seconds after GW170817 
2017). The occurrence of an electromagnetic counterpart — the kilonova — is a 
great advantage of a DNS merger or a NS+BH merger, relative to the mergers 
of double BHs, because this allows one to determine the place of origin of the 
merger on the sky with sub-arcsecond precision. In the case of GW170817, this 
place a EE out to 2 in the lenticular (S0) galaxy, NGC 4993 at a distance of 
40 Mpc (C [Coulter et al.|[2017} a (2017). The discovery, in this 
way, of the DNS merger GW170817 at once provided the solution to the nature 
of a variety of key astrophysical phenomena. 

In the past decades it was realized that the SNe of Types Ia, Ib and Ic 
(characterized by the absence of hydrogen in their spectra) most likely are re- 
lated to the evolution of binary systems (see Section 2.10). That these three 
types together form about half of all SNe shows how important knowledge of 
binary star evolution is for a general understanding of observational stellar evo- 
lution processes. 

Finally, in the past decades it was discovered that binary evolution has 
affected a large variety of low- and intermediate-mass stars with peculiar abun- 
dances of chemical elements. Examples are the barium stars, which are G and 
K giants with masses up to a few solar masses and an overabundance of bar- 
ium and other s-process elements (products of slow neutron-capture processes 


(Bidelman & Keenan} |1951 Bidelman & Keenan} {1951} [Boffin & Jorissen] Boffin & Jorissen\ |1988 [1988). They are wide binaries 


with eccentric orbits and orbital periods ranging from approximately 100 to 
10000 days (1988). The binary nature of these stars, of blue 
stragglers in star clusters, of extremely metal-depleted post AGB stars, of car- 
bon enhanced metal poor (CEMP) stars, and of over 80% of nuclei of planetary 
nebulae was discovered in the past decades. These discoveries show that the 
evolution of ~50% of all low-and intermediate-mass stars is expected to be af- 
fected by binary evolution 2015). Because lower-mass stars are 
far more abundant in galaxies than massive stars are, the products of low- and 
intermediate-mass binaries are expected to vastly outnumber the products of 


the evolution of massive binaries. However, in this book, we will concentrate 
mainly on the physics and evolution of the more massive binaries that produce 
NSs and BHs and on the origin of Type Ia SNe, which are caused by exploding 
WDs and are of crucial importance for cosmology. 

In the last decades it was also realized that stars are found often even in 
triple systems or higher-order multiple systems. Observational estimates sug- 
gest that approximately 20 — 30% of all binary stars are members of triple 
systems BOI). These systems may 
remain bound with a long-term stability in a hierarchical structure (a close inner 
binary with a third star in relatively distant orbit). Since the 1970s, a number 
of stability criteria for a triple system have been proposed (see 2008) 
for an overview) that enable predictions for their long-term stability. 
estimated that in ~70% of the triple systems, the inner binary 
is close enough that the most massive star will evolve to fill its Roche lobe. 
Furthermore, in ~ 15% of the triple star systems, the outer third (tertiary) star 
may even also fill its Roche lobe at some point, possibly leading to disintegration 
or production of rare configurations with three degenerate objects in the same 
system. In 2014, reported the remarkable discovery of a 
triple-system pulsar (PSR J0337+1715), which is exactly the first example of 
such an exotic system — a NS orbited by two WDs. (For a formation and evolu- 
tion scenario for this complex system, which must have survived a SN explosion 
and at least three stages of mass transfer between the stellar components, see 
our model in [Tauris & van den Heuvel] [2 2014). Besides specific ual star in- 
teractions, such as Kozai-Lidov resonances (Kozail |1962 Kozai] [1962] [Lidov} [1962] 1962) and the 
above-mentioned dynamical stability considerations, most a aon between 
stars in triples are similar to the interactions between binary stars. For this 
reason, we focus only on binary stars in this book? | 

The 2020s and 2030s are expected to reveal a large number of discoveries 
of new double compact object systems, as well as their progenitors and merger 
remnants. This field of astrophysics is strongly driven by investments in new big- 
science instruments. The Square-Kilometer Array (SKA) is expected to increase 
the number of known radio pulsars o a T of 5 to 10, thus resulting in a to- 
tal of > 100 known DNS systems Keane et al.|{2015). The Five-hundred-meter 
Aperture Spherical Telescope SHa is alo ETE to a. a significant 
number of new radio pulsars (S (Smits et al.|{2009}[Nan et al.}[2011) , including new 
discoveries of pulsar binaries. High-mass X-ray binaries (HMXBs), the antici- 
pated progenitors of double compact object systems containing NSs and BHs, 
are continuously being discovered with ongoing X-ray missions (INTEGRAL, 
Swift, XMM-Newton, and Chandra, see e.g. (2013). New and upcom- 
ing space-borne X-ray telescopes such as eX TP, STROBE-X, and Athena are 
expected to produce further discoveries of these systems. Hence, we are cur- 
rently in an epoch when a large wealth of new information on exotic binaries is 
becoming available. In light of this, it is important to explore and understand 


?Formation of stellar binaries via triple- or quadruple-star dynamical interactions in dense 
cluster, however, are discussed in Chapter 12. For a recent review on the evolution of desta- 


bilized triple systems, see{Toonen et al.| (2022). 


the formation and evolution of such binary systems in more detail. Earlier 


textbooks on binary evolution are those of (1994), and |Egeleton 


(2006), to which we refer for further reading. For an earlier book on physics 
and evolution of relativistic objects in binaries, we refer to [Colpi et al.] (2009). 

This book is organized as follows. In Chapter 2, we give a brief history of 
the discovery of the different types of binary systems and, where appropriate, 
summarize their importance in modern astrophysics. In Chapter 3, we consider 
how the orbital parameters of spectroscopic and eclipsing binaries are measured 
and how from these measurements, information is obtained about the masses 
and radii of the stars. We give an overview of the thus derived masses of stars 
of different spectral types. 

In Chapter 4, we consider basic aspects of the celestial mechanics of binary 
systems, the meaning and limitations of the Roche-lobe concept, as well as the 
changes in orbital period and binary separation that are induced by various 
processes of mass loss and mass transfer in binary systems. We first consider 
the somewhat idealized “conservative” evolution, in which the total mass and 
orbital angular momentum of the binary are expected to be conserved during 
the evolution, followed by the more realistic and exotic types of close binary 
evolution, so-called non-conservative evolution, in which large losses of mass 
and orbital angular momentum from the systems are taken into account. A 
treatment of common envelopes and the orbital evolution during the dynami- 
cally unstable in-spiral phase are described as well. Finally, we briefly discuss 
the Eddington accretion limit as well as accretion disks. 

In Chapter 5, we describe the observed properties and the general classifi- 
cations of the various types of interacting binary systems, that do not contain 
NSs or BHs, concentrating mostly on systems in which at least one component 
is an evolved star, that is, a (sub)giant, or a WD. 

In Chapters 6 and 7, we describe the observed properties of X-ray binaries: 
high-mass X-ray binaries (HMXBs) as well as low-mass X-ray binaries (LMXBs), 
including mass determination of the accreting NSs and BHs. Regarding HMXBs, 
we discuss Be-star X-ray binaries, supergiant X-ray binaries, and, for example, 
stellar wind accretion and the Corbet diagram. We also discuss the recently 
discovered class of pulsating ultra-luminous X-ray sources. For LMXBs, we 
discuss the various types, including the systems with BHs and the symbiotic 
X-ray binaries with accreting NSs. 

In Chapter 8, we give an overview of the evolution of single stars (with a 
special focus on the final evolution of massive stars). In Chapter 9, we apply this 
knowledge to the evolution of binaries in general. Here, we also discuss in detail 
the various cases of mass transfer (including mass loss) and orbital stability 
analysis, and we end with a comparison between the outcomes of single versus 
binary star evolution. 

This knowledge is crucial for understanding the formation and evolution 
of X-ray binaries, which are the subjects of Chapters 10 (HMXBs) and 11 
(LMXBs). We discuss the final stages of HMXBs (including Wolf-Rayet star 
binaries), with or without a common envelope, leading to the formation of dou- 
ble NS/BH binaries. For LMXBs and CVs, we also discuss the mechanisms 


driving the mass transfer in LMXBs and CVs. These concern the internal evo- 
lution of the companion star, as well as the loss of orbital angular momentum 
due to emission of GWs and/or a magnetically coupled stellar wind, or a combi- 
nation of these. We also discuss the final mass-transfer stage from WDs in very 
tight binaries, in the so-called AM Canum Venaticorum (double WD) systems 
and the ultra-compact X-ray binary sources (typically WD+NS systems). 

Binaries with compact objects can also be formed by the dynamical evolution 
of dense star clusters. This is the subject of Chapter 12. 

Chapter 13 concerns SNe in binaries. We first discuss the evolution leading 
to the thermonuclear Type Ia SNe, triggered by the accretion of matter by 
a WD, and we also consider accretion-induced collapse of massive WDs. We 
subsequently discuss the origin of Type Ib/Ic SNe, and the evidence derived from 
theoretical computations of the late stages of close binary evolution, including 
ultra-stripped SNe, combined with observations of SN light curves and known 
Galactic post-SN DNS systems. We discuss the evidence for momentum kicks 
imparted onto compact objects during various SN events (iron core collapse and 
electron-capture SNe) and examine the kinematical effects of these kicks on the 
resulting compact-object binaries. 

In Chapter 14, using the results of the previous chapters, the final stages of 
the formation of binary and millisecond radio pulsars is studied; in other words, 
the transition from X-ray binaries to radio pulsars. We review the recycling of 
pulsars in detail, including the accretion torques at work. We review the rich 
diversity of resulting millisecond pulsar binaries and their component masses 
(theoretical expectations and measurements via relativistic effects). Finally, we 
revisit the formation of DNS systems. 

Chapter 15 is devoted to GW astrophysics. We cover the basic physics of GW 
emission and detection. We discuss formation channels of both high-frequency 
(LIGO-Virgo-KAGRA-IndIGO) and low-frequency (LISA and TianQin) GW 
sources. We review the signals expected from extragalactic merging double BHs 
and DNSs, including the electromagnetic counterparts of these mergers, and the 
Galactic WD and NS binaries as continuous GW emission sources. The different 
models for the formation of the double BHs are discussed in depth in light of 
the latest observations of “ordinary” and exotic events from the LIGO-Virgo-— 
KAGRA network. 

Chapter 16 discusses the subject of binary population synthesis with an 
emphasis on methodology and statistics. Two examples are highlighted that 
illustrate a synthetic open star cluster population of binaries and the differences 
between estimates of empirical versus theoretical DNS merger rates. 
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